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1. Introduction 



Many new states have been observed in the charmonium sector in the past few years. The miss- 
ing spin singlet states as well as new, yet to be explained states such as the X(3872)[jl|], Y(4260)[§], 
Y(3940)[j3|] and others have presented a new challenge for the theoretical particle physics commu- 
nity. 

Simulating the charmonium spectrum in lattice QCD has long been problematic due to the 
fact that the mass is too heavy for a simple relativistic action on coarse lattices and too light to 
expand in inverse powers of the quark mass. There are various quenched studies with improved 
actions which attempt to reduce lattice artifacts [Q, ||] but very little has been done on dynamical 
configurations 

In this study, we use a new dynamical, anisotropic action which has the advantage of being 
fully relativistic and yet has small mass-dependent discretization errors. All-to-all propagators [^] 
are utilized to make full use of the configurations and to construct extended operators for better 
overlap with the higher charmonium states. 



2. Dynamical anisotropic action 

The Nf = 2 anisotropic quark action is as follows: 



Sq = ¥ f %V + £ HryVi (l - |~2 A ') - "jjf A '0 + L A ^ + m o^J V , ( 2 - 

where the links are fattened using stout links [^] which maximize the plaquette. The target renor- 
malized anisotropy, £ , is 6. The nonperturbative tuning of the action has been presented elsewhere 
at this conference []10||. The sea quark mass was set near the strange quark in this first study al- 
though there are indications that there are no major problems simulating at lighter sea quark masses. 



The gauge action is the "two-plaquette" action [ |110 which was designed to reduce the cutoff effects 



for the scalar glueball on coarse anisotropic lattices. We use the same anisotropic fermion action 
for both light sea quarks and heavy valence quarks. In a quenched study it was found that the 
same bare anisotropy can be used for a large range of quark masses. For the heavy valence quark, 
we have used both stout links and original links in this study. 



3. Operators 

We use a variational basis of operators, many of which are extended in space, to maximize 
the overlap with the higher charmonium states. All-to-all propagators which do not introduce 
large variance are crucial for this construction. We use the "dilution" method of Ref. [Q] without 
eigenvectors for the charm quark propagators. Several dilution schemes were tested on a small 
number of configurations. We present results which use "time+colour+space-even-odd" dilution. 
Preliminary investigations indicate that even higher dilutions may further reduce the errors; this is 
work in progress. 

The S-wave charmonium states only require different smearings of the quark to project out 
radially excited states. However, we have included a hybrid 1 operator which contains a gluonic 
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Particle Operator 



1 f, 75«> YsB 

0~ + 75,75(^1+^2 + ^3) 

0++ y.p 

1 ++ Y X P 

2 ++ JkPi + 7>Pk 

l + ~ 75^ 

1~ + yxw, yxfi 

2~ + 75(^1 -53), 75( 2 ^3 -S3) 

3~~ y-t 



Table 1: Some operators used to project out various higher lying states. The notation for the gluonic paths 
isthatofRef. @. 



excitation in the form of a chromomagnetic and "staple" field. The operators used in this study 
are summarized in Table [I]. Many of the operators used here resemble the ones used in Ref. [jTJ]. 
In addition to these basic operators, we have used two different smearings for the quark fields to 
build a larger variational basis. The variational optimization was performed at the largest timeslice 
possible (with the metric timeslice fixed at t = 1) where the results have become independent of 
this choice. We note that the r] c and % c states (especially the ++ ) required a large t for the results 
to become stable. This may be an indication that better operators can improve the situation for 
these channels significantly. 

For the exotic channel, 1~ + , we have used both the staple and chromomagnetic field operator 
since we know that the staple operator has a large overlap with the first excited gluonic state in 
the static qq system [13]. The chromomagnetic field operator has the advantage of being an easy 
operator to construct as it is an operator at a single point. We have found that there is very little 
difference between the two operators, contrary to the non-relativistic quark model which suggests 
that the quarks are in a relative P-wave. This may be due to the fact that the quarks have been 
smeared or that both operators have a high level of excited state contamination from multi-particle 
states. 



4. Analysis 

4.1 Fitting scheme 

The use of time-diluted all-to-all propagators makes it difficult to identify the plateau region 
(and judge its quality) from effective mass plots. This is due to the random noise introduced at each 
timeslice which makes locally measured quantities such as the effective mass fluctuate more than 
they would with point propagators. These "local fluctuations" do not affect exponential fits as long 
as a reasonable range of timeslices are included since these fits capture the long range exponential 
decay of the correlation function. (This feature will disappear once the gauge noise limit is reached 
in the dilution method. We already see a noticeable decrease of these effects at the higher dilution 
results that are presented here.) We therefore show instead "? m i n " plots where fitted values of the 



029/3 



Charmonium Spectrum 



K. Jimmy Juge 



Splitting 


stout 


mo 




ai l (GeV) 


P-S 


with 


0.09 


0.0568(5) 


8.06(7) 


P-S 


w/o 


0.12 


0.0668(7) 


6.85(7) 



Table 2: Temporal lattice spacings determined from spin-averaged SP splittings in the charmonium sector 
for the two actions used in this study. 

mass are shown against the minimum timeslice included in the single exponential fits. The largest 
timeslice to include in the fit was taken to be halfway across the lattice, if the fluctuations were not 
dominating the signal. We show a typical example of an effective mass (1 ++ ) and its corresponding 
tmm plot on the left. Our criteria for selecting a fit is based on the # 2 , fit range and the fit quality 
Q. We demand (although there could have been many other choices for the criteria for a "good fit") 
that X ,/Ndf be less than 2, the quality of the fit be larger than 0.2 and that there is a range of t m \ n 
where these fit values are stable. 
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Figure 1: The effective mass of the 1 ++ ground state is shown in the left figure. The "tmm plot" for this state 
is shown on the right. The open symbols in this plot indicate that the fit did not pass the requirements for a 
good fit (see text). 



4.2 Setting the scale 

We estimate the temporal lattice spacing from the S-P splitting in the charmonium sector. The 
results are tabulated in Table ||. We note that the bare quark masses were not quite tuned to the 
charm quark mass and the puzzling dependence on the mass and/or the effects of the stout-link is 
still under investigation. 

4.3 Dispersion relation 

In order to check that the aM discretization errors are under control, we compute the dispersion 



relation for r\ c (prJ/xjf) at various values of the momenta. The result is shown in Fig. 4.3. Since the 
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renormalized anisotropy enters into this calculation, we cannot at this point determine unambigu- 
ously whether the kinetic mass is equal to the rest mass or not. This would require a determination 
of the dispersion relation of another state, e.g. D s . We find a renormalized anisotropy determined 
at the charm quark mass which is significantly higher than what is found at the strange quark mass 
[14]. This discrepancy may be due to aM errors, but it may also have its origins in the different 
methods used to determine the anisotropy (point propagators in [|l4|] vs. all-to-all propagators here), 
combined with the short temporal extent of our lattice. This is currently under investigation, using 
lattices with a longer time extent. 




5. Results and discussion 

The spectrum computed with (without) stout links is shown on the left (right) of Fig. ||. Since 
the quark mass has not been exactly tuned to the charm mass, the spectra are shown with the mass 
of the J/\\f adjusted to its physical value. The main difference between the two spectra appears 
to be a change in the scale. Whether this is an artifact of the fitting criteria or the diagonalization 
procedure, or whether it is a real effect, is under investigation. 

5.1 Hyperfine splitting 

All previous lattice studies have given a too small value for the hyperfine splitting. It is clear 
that this problem is not yet resolved with this action. At this stage, the T] c and J/y are degenerate 
within the errors. There are several reasons why this splitting is not observed: 

• Finite volume: The sides of this lattice are roughly 1.6 fm, which is quite small for charmo- 



nium. At this volume, finite temperature effects may also play a role [jlOj. 



Higher order improvements of the action: The action used in this study has leading dis- 
cretization errors of ff{a t ,a s a S: cr s ). Higher order correction terms, in particular the a B 
term, may have a significant impact on the hyperfine splitting. 
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Figure 3: The charmonium spectrum from the simulation using stout links (left) and without stout links 
(right). The quantum numbers of some of the experimental findings for the higher lying P-waves have not 
been established yet; the assignment shown here is the one presented in Ref. Jl5[], 

• Light sea quark effects: The fact that the hyperfine splitting comes out too small in the 
quenched approximation indicates that light sea quark effects are important, and the sea 
quark mass that we are using may still be too large. 

5.2 States above threshold 

We have computed the masses of states which may lie above threshold, even at the relatively 
heavy sea quark mass used here. If this is so, then we must include multi-particle states in the 
spectrum and perform a finite volume analysis before we can claim that we have determined the 
energy of these states. Multi-particle operators are accessible with all-to-all quark propagators and 
investigations in this direction are under way. 

Figure [| shows the effective mass for the exotic hybrid (1~ + ), which is found to lie higher 



than in some quenched calculations (see [ |16| ] for example). For this case, however, it is probably 
even more important that a larger volume is used since the wavefunction of the exotic hybrid could 
be quite large [ |T7| ] and would be squeezed in the small lattices used here. 

5.3 Hybrid content of 1 

The correlation matrix element in the 1 channel between the state created using a standard 
eye operator and a hybrid cy^Bc operator is plotted in Fig. ||. The coefficient appears to be zero 
within errors, suggesting that the hybrid content of the y spectrum is small. The hybrid content 



of the T spectrum was also found to be small in a quenched calculation [18]. However, further 
reduction in noise is required before drawing any definitive conclusion. 

5.4 Disconnected diagrams 

Disconnected diagrams were not evaluated in this first study as these contributions are ex- 



pected to be small (see [19] for example). However, the inclusion of these diagrams is not difficult 
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Figure 4: Effective mass of the 1 
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Figure 5: Correlation functions of the 1 state with different interpolating operators. 



with the all-to-all propagators used here. A complete study with tuned parameters is planned to 
study the effect of these neglected diagrams on hyperfine splittings. 

6. Outlook 

A first study of the charmonium spectrum from an anisotropic, Nf = 2 action was presented. 
All-to-all propagators and many operators were used to optimize the overlap with the standard 
states and obtain signals for higher lying radial excitations. A good signal for D-waves and exotic 
hybrid states was obtained with a relatively low level of dilution. However, short distance effects 
such as the hyperfine splittings require further investigation. A more detailed study on a larger 
volume and with more precisely tuned parameters is in progress. 
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